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Ph. On the basis of the results of a new renormalisation group improved small- x resumma- 

<-j ! tion scheme, we argue that the range of validity of perturbative calculations is considerably 

extended in rapidity with respect to leading log expectations. We thus provide predic- 
tions for the energy dependence of the gluon Green function in its perturbative domain 
and for the resummed splitting function. As in previous analyses, high-energy exponents 
are reduced to phenomenologically acceptable values. Additionally, interesting preasymp- 
totic effects are observed. In particular, the splitting function shows a shallow dip in the 
moderate small- x region, followed by the expected power increase. 
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1 Introduction 



The prediction of high-energy hard cross sections in QCD perturbation theory has been a 
puzzling problem in recent years for a number of reasons: the existence of large perturbative 
leading logs (LL) contributions pQ that seem incompatible with a range of experimental data; 
the discovery of even larger subleading contributions (NLL) of opposite sign pHE]; and the 
increasing importance of low-/c t partons at high energies leading asymptotically to a strong- 
coupling Pomeron regime which can at best be modelled but not really predicted jHEHEllZl • 

In order to tame the problem of large logarithms with alternating sign several resummation 
procedures have been proposed [3131101111111211131111113 • The renormalisation group improved 
(RGI) approach |H|PITU] offers a general understanding of the large subleading log s terms as due 
to leading log Q 2 collinear contributions which can be extracted on the basis of the perturbative 
anomalous dimensions. As a result, resummed high-energy exponents have been calculated in 
a stable way In a forthcoming paper PB] this resummation approach is extended to the gluon 
Green's function and splitting function so as to provide a fairly complete account of their energy 
dependence and of its perturbative (PT) versus non-perturbative (NP) aspects. The purpose 
of the present paper is to summarise the main results of Ref. [TBI and to argue, on this basis, 
that the RGI approach tames - to a large extent - the problem of the strong-coupling region 
as well. This is because the resummation of the (alternating sign) large logarithms leads to 
smaller high-energy exponents (and diffusion coefficients) and to a considerable suppression of 
the non-perturbative Pomeron itself in the Green function. Furthermore, we find that strong- 
coupling contributions factorise, as expected, in the collinear limit, and we are able to provide 
the resummed perturbative splitting function in x-space. 

The basic problem that we consider is the calculation of the (azimuthally averaged) gluon 
Green function G(Y; k, ko) as a function of the magnitudes of the external gluon transverse 
momenta k = \k\, k = \k \ and of the rapidity Y = log -g^. This is not yet a hard cross section, 
because we need to incorporate [IB] the impact factors of the probes [T71ITH1IIII] • Nevertheless, 
the Green function exhibits most of the physical features of the hard process, if we think of 
k 2 , as external (hard) scales. The limits k 2 3> k^ (k^ 3> k 2 ) correspond conventionally to 
the ordered (anti-ordered) collinear limit. By definition, in the tu-space conjugate to Y (so that 
Cj = dy) we set 

g u (k,ko) = [w-JCJ-^fco), (1) 
u>g u (k, k ) = 5 2 (k -ko) + J d 2 k' K u {k, k')Qu{k', k ) , (2) 

where /C w (fc, k') is a kernel to be defined, whose u = limit is related to the BFKL ^-evolution 
kernel, known at LL and NLL levels [21E1- 

The RGI approach is based on the simple observation that, in BFKL iteration, all possible 
orderings of transverse momenta are to be included, the ordered (anti-ordered) sequence k ^> 
ki ■ • • ^> k n • • • ^> k Q {k <C k\ • • • k n • ■ ■ k ) showing scaling violations with Bjorken 
variable k 2 / s (k^/s). Therefore, if only leading log A; 2 contributions were to be considered, the 

kernel acting on p- ( J would be approximately represented by the following eigenvalue 



2 



function (in the frozen coupling limit) 



o:. s 



10 



7 + 



1 + 



7 



10 



+ A 1 (oo) + 



Nr. 



O.S 



a. 



71 



(3) 



where 7^ = a s (~ + Ai(u)) is the one-loop gluon-gluon anomalous dimension and we have 
introduced the variable 7 conjugate to logfc 2 . Note, in fact, that Eq. (jSJ) reduces to the normal 
DGLAP evolution [20] in log/c 2 (log A; 2 ,) in the two orderings mentioned before, because 7 + ^ 
(1 + ^ — 7) is represented by <9 logfc 2 (<9 logfc 2) at fixed values of x — k 2 /s (x = k^/s) in the 
ordered (anti-ordered) momentum region. Note also the tu-dependent shift (H1EIJIIII! of the 7- 
singularities occurring in Eq. (jnj), which is required by the change of scale (fcg versus k 2 ) needed 
to interchange the orderings, i.e., x versus x. 

We thus understand that the tu-dependence of fC^ is essential for the resummation of the 
collinear terms and can be used to incorporate the exact LL collinear behaviour, while on the 
other hand, the uj —>■ behaviour of /C w is fixed by the BFKL limit up to O (uj) terms, so as to 
incorporate exact LL and NLL kernels. Such requirements fix the kernel up to contributions 
that are NNL in In a; and NL in InQ 2 . 1 The form proposed in ^B] (there called NLLb) is given 
by 



= {a s (xlq 2 ) + ba 2 s In zj) Kq + u (a s (xjA; 2 ) + ba\ In x 2 ) K? + a 2 ^ 2 )K? 



(4) 



x M has been introduced to test 



where q = \k — k'\, fc> = max(fc, k'), k< = mm(k,k') 
renormalisation-scale dependence, and the scale-invariant kernels in the RHS have the eigen- 
value functions 



X£( 7 ) = 2^(1) - ^(7 + |)-V(l-7+f), 
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with b = li - C(u;) = ili±^ffl _ an d 

12 6iV c ' \ / oi a;+l 



Xo(7) = Xo = °(7) = 2^(1) - V(7) - ^(1 " 7) 



Xi(7) = -n[Xo(7) + Xo(T)] - 7X0(7) - 7 
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1 We choose not to include the known exact NL InQ 2 terms because the conceptual issues, associated for 
example with factorisation scheme dependence, have yet to be fully understood. 
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All rif dependence other than that in the running coupling terms (that is in b) has been 
neglected. A correct treatment of all nj-dependent RG contributions would make the formalism 
technically more complex (e.g. requiring a two-channel approach) and given the observed small 
effect of the n/ contributions, we feel such an effort to be currently uncalled for. 

Note that, because of the explicit form of Xo , the kernel (jH) reproduces the pole behaviour (J3J) 
to first order in a s . Note also that the last line in Eq. (JScJ), which vanishes in the u = limit, 
has been added in order to shift the remaining C(0)/j and C(0)/(1 — 7) poles in \i- There is 
some freedom in the choice of the a;- dependence of the coefficient of the shifted poles, and we 
take here a minimal prescription, one that gives an indentically zero two-loop contribution to 
the anomalous dimension (we refer to this prescription as NLLb). This also guarantees that 
the momentum sum-rule is satisfied at two-loop level. Note finally that the resummed kernel 
proposed here differs from that of [[TO] because the collinear terms are added in the a;-dependent 
form, the difference being at NNLL level. The reason for such a change is that we have at most 
simple collinear poles, as in Eq. (|SJ), and not double poles. This avoids the need for the u> 
expansion, thus providing a kernel more suitable for numerical iteration in (x, fc)-space. 

The integral equation to be solved by the definition © is thus a running coupling equation 
with non linear dependence on a s at appropriate scales, and it has a somewhat involved in- 
dependence in the improved kernels (J5J). Its solution has been found in JH] by numerical matrix 
evolution methods in k- and x- space |21j . where the typical u;-shifted form in the example Q 
corresponds to the so-called consistency constraint [22 ESI El] • Furthermore, introducing the 
integrated gluon density 



and has been extracted ^Hj by a numerical deconvolution method 25j. We note immediately 
that P e s turns out to be independent of k for Q 2 ^> k 2 ,, yielding an important check of RG 
factorisation in our approach. 

2 Gluon Green's function 

Results for G(Y; k, k) 2 are shown in Fig. ^ In addition to the solution based on our RGI 
approach (NLLb) the figure also has 'reference' results for LL evolution with kernel a^x^q 2 )^. 
We use a one-loop running coupling with nj = 4, which is regularised either by setting it to zero 
below a scale k ('cutoff') or by freezing it below that scale (a s (q 2 < k 2 ) = ot s (k 2 )). We believe 
the cutoff regularisation to be more physical since it prevents diffusion to arbitrarily small scales 

2 Actually, we take slightly different values of the scales in G, namely we consider G(Y; k + e, k — e) with 
e = O.lfc, in order to avoid sensitivity to the discretisation of the (5-function initial condition in J2J) (cf. Ref. ^B] 
for a detailed discussion). 





the resummed splitting function P e s(z,Q 2 ) is defined by the evolution equation 




(9) 
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Figure 1: Green's function calculated with four different infrared regularisations of the coupling, 
shown for LL and RGI NLL ('NLLb') evolution. The bands indicate the sensitivity of the 
k = 0.74 GeV results to a variation of x 2 ^ in the range | to 2. The left and right hand plots 
differ only in their scales. 



and is thus more consistent with confinement — accordingly we show three cutoff regularisations 
and only one frozen regularisation. The k = 0.74 GeV cutoff solution is presented together with 
an uncertainty band associated with the variation of x 2 ^ between | and 2. 

Solutions of (0) with an IR- regularised coupling generally have two domains |H 03 01211], 
separated by a critical rapidity Y c (k 2 ). For the intermediate high-energy region 1 < V < Y c (k 2 ), 
one expects the perturbative 'hard Pomeron' behaviour with exponent u s , 

k 2 G(Y; k,k)~^= exp [u s (a s (k 2 ))Y + A(a s , Y)] , (10) 
V Y 

and diffusion corrections |2T | l2l ^ l2T H I!?n] parametrised by A(a s , Y). Beyond Y c , a regularisation- 
dependent non-perturbative 'Pomeron' regime takes over 

J I.., (ID 

where the factor £ differs from 1 only for kernels involving the consistency constraint. The 
non-perturbative exponent Uj? satisfies [TU] cup ~ u s (a s (k 2 )) and hence is formally larger than 
^,(n,.(/. /J )!. :i 

The value of Y c depends strongly on k. In the tunnelling approximation, it can be roughly 
estimated by equating eqs. (fTU|) and (JTTJ) to yield [HUEE], for any given regularisation procedure, 

3 The behaviour (|llfl with u>p > uj s (a s (k 2 )) is a general feature of linear evolution equations such as (0), but 
not of actual high energy cross sections, which are additionally subject to non-linear effects and confinement. 
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u F - u s (a s (k 2 )) 

(again with 1 + up — > 1 for LL), showing an approximately linear increase of Y c with logfc 2 . 

Within this logic, several aspects of Fig. are worth commenting. The most striking feature 
of the LL evolution is its strong dependence on the non-perturbative regularisation, even for 
rapidities as low as 5. The exact value of Y c depends on the regularisation being used, ranging 
between 5 and 10. In contrast, NLLb evolution remains under perturbative control up to much 
larger rapidities and the NP pomeron behaviour takes over only for Y > 25, where the three 
cutoff solutions start to diverge. Regularisation dependence is also present at lower rapidities. 
This seems to be due to power corrections to (fTUJ) . associated with the use of a coupling 
a s(q 2 ) P^l - The non-perturbative regime of the IR frozen coupling solution is reached only 
later (Y > 30), at the point where it starts to grow more rapidly than the cutoff solutions. A 
final point to note is renormalisation scale uncertainties of our resummed results are sizeable - 
of the order of several tens of percent for Y > 4 - but seem anyhow quite modest compared to 
the order(s) of magnitude difference with LL. 4 

The fact that Y c is considerably larger for NLLb evolution is natural — it is an expected 
consequence (see Eq.(|12Jl) of the fact that subleading corrections lower both the PT and NP 
exponents. What is quite non-trivial however is that at their respective Y" c 's the NLLb Green's 
function is an order of magnitude larger than the LL one: the subleading corrections increase 
the overall amount of BFKL growth remaining within perturbative control. This is due to a 
number of factors, among them a strong reduction of the diffusion coefficient (see below). We 
note that the increase in the amount of perturbatively calculable growth is of particular interest 
for the theoretical question (see e.g. |2Zj) of whether it is possible to perturbatively generate a 
high-density gluonic system. 

In Fig. ^ we have considered only a single value of k. The question of NP contributions 
is summarised more generally in Fig. which shows contour plots of the logarithmic spread 
of our four regularisations. Darker regions are less IR sensitive, and contours for particular 
values of the spread have been added to guide the eye. Here too one clearly sees the much 
larger region (including most of the phenomenologically interesting domain) that is accessible 
perturbatively after accounting for subleading corrections. 

So let us now therefore return to Fig.^and examine the characteristics of the NLLb Green's 
function in the perturbatively accessible domain, which should be describable by an equation 
of the form (|10|). The first feature to note is that the growth starts only from Y ~ 4. This 
suggests that at today's collider energies (implying Y < 6 [3^113^] ). it will at best be possible to 
see only the start of any growth. This preasymptotic feature is partly due to the slow opening 
of small-x phase space [34] implicit in our cu-shifting procedure. 

Once the growth sets in, the issue is to establish the value of u> s appearing in (fTUjl. This 
is a conceptually complex question because in contrast to the fixed-coupling case, uj s no 
longer corresponds to a Regge singularity. There are running-coupling diffusion corrections 

4 The x^ dependence of the LL solution is somewhat smaller than for NLLb — this may seem surprising, but 
at larger Y the LL solution is in the NP domain, where non-linearities (in a s ) reduce the x^ dependence. 
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Figure 2: Contour plots showing the sensitivity of G(Y, k + e, k — e) to one's choice of non- 
perturbative regularisation, as obtained by examining the logarithm of the ratio of the regulari- 
sations giving the largest and smallest result for G. Darker shades indicate insensitivity to the 
NP regularisation, and contours have been drawn where the logarithm of the ratio is equal to 
0.1, 0.2 and OA. Plot (a) shows the result for LL evolution, while (b) shows RGI NLL evolution 
(NLLb). The regularisations considered are those of Fig. 



A(a s ,Y) |2Zll2Hll2nilSni , whose leading contribution, ~ Y 3 , for our model is [THj 

as) 

In addition, A(a s ,Y) contains terms with weaker Y dependences, including Y 2 and Y. Such 
terms can be disentangled by the method of the 6-expansion [35 . 5 Since running-coupling 
diffusion corrections start only at order b 2 , it is possible to give an unambiguous definition of 
uj s up to first order in b, while retaining all orders in a s for non running-coupling effects. The 
result is shown in Fig. El as a function of a s and as has been found in previous work ^U|, there 
is a sizeable decrease with respect to LL expectations. 6 Furthermore, the leading diffusion 
corrections in (jlHj) turn out to be numerically small, about an order of magnitude down with 
respect to the LL result, due to a sizeable decrease of the diffusion coefficient x'^si over an d 
above the decrease already discussed for u s . 

A final point related to the Green's function concerns predictions based on a pure NLL kernel 

5 In the d -t limit, with a s (k 2 ) kept almost fixed, the non-perturbative Pomeron is exponentially sup- 
pressed |H5H12j . so that the 6-expansion can also be used as a way of defining a purely perturbative Green's 
function without recourse to any particular infrared regularisation of the coupling . 

6 We do not compare directly with our earlier results for cu s I10j . because they are based on a different 
definition (the saddle-point of an effective characteristic function), which is less directly related to the Green's 
function. Nevertheless, the present results are consistent with previous ones to within NNLL uncertainties. 
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Figure 3: The small-x exponents: the Green's function effective exponent u s is shown to first 
order in the b- expansion; the splitting function exponent u c is shown together with NP and 
renormalisation scale uncertainty bands, defined in figure^ Also shown, for reference, is the 
result for u c using the method of ITf^ . for 6(n/ = 4). 

without renormalisation group improvement. Despite the fact that the characteristic function 
around 7 = ~ is very poorly convergent, it has been argued [SEj that the Green's function may 
nevertheless show a growth corresponding to an effective positive u s , due to saddle-points at 
complex 7. Indeed, we find jTH] that for k ~ ko the pure NLL Green's function is remarkably 
similar to the NLLb result and it is stable with respect to x^ variations. There is a difference 
in normalisation, which we ascribe to to the (implicit) presence of an effective impact factor in 
the NLLb solution [T^]. It should be kept in mind however that this 'good behaviour' of the 
pure NLL Green's function breaks down when k and ko are substantially different (for a s ~ 0.2 
the Green's function becomes negative for | \nk/ko\ > 2), resulting in the unphysical oscillating 
behaviour predicted in [36J. 

3 Resummed gluon splitting function 

Green's functions G(Y; k, ko) with k ~ k ~ Aqcd are very sensitive to, and largely determined 
by non-perturbative physics associated, in our numerical solutions, with an IR regularisation 
of the running coupling at some scale k. For example at x = e~ Y = 10~ 10 and k = 4.5 GeV, 
the three cutoff regularisations of the previous section lead to a spread of a factor of 5 in 
calculations of the Green's function, in sharp contrast to the good perturbative control seen in 
the previous section, for the case ko ~ k ^> Aqcd- 

However, many arguments in the BFKL framework, jSJHJITJEIlEllSIiL have been given in 
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favour of factorisation, Eq. Q, with the small-x splitting function P gg (z, Q 2 ) being independent 
of the IR regularisation. The most dramatic demonstration of factorisation is perhaps in the fact 
that a numerical extraction of the splitting function from the Green's function by deconvolution 
gives almost identical splitting functions regardless of the regularisation. This is illustrated in 
Fig. HJ where the solid line and its inner band represent the result of the deconvolution together 
with the uncertainty resulting from the differences between the three cutoff regularisations. 7 
The resulting regularisation dependence is pretty small, and at higher Q it diminishes rapidly 
as an inverse power of Q, as expected from a higher twist effect. 
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Figure 4: Small-x NLLb resummed splitting function, compared to the pure 1-loop DGLAP 
and the (fixed- coupling) LL BFKL splitting functions. The central NLLb result corresponds to 
Xfj, — 1, k — 0.74 GeV; the inner band is that obtained by varying k between 0.5 and 1.0 GeV, 
while the outer band corresponds to \ < x 2 ^ < 2. 

Several features of the resummed splitting function are worth commenting and comparing 
with previous NLL calculations, including various types of resummations [lT H IT2 t lT8] ■ Firstly, 
at very large-a; it approaches the normal DGLAP splitting function. The momentum sum rule 
is satisfied to within a few parts in 10 4 . At moderately small-x the splitting function is quite 
strongly suppressed with respect to the LL result and shows not a power growth, but instead a 
significant dip (of about 30% relative to the LO DGLAP value, for a s = 0.215). Dips of various 
sizes and positions have been observed before in [T8|llip25j though ours is significantly shallower 
than that found in ^3] at NLL order and similar to that found in the cu-expansion ^Bj- An 
interesting question concerns the impact of the dip on fits to parton distributions. Calculations 
in a (partially) RGI LL model [21] whose effective splitting function also has a dip, suggest 

Because of numerical instabilities in the rather delicate deconvolution procedure, we have so far not suc- 
ceeded in obtaining reliable results with an IR frozen coupling — preliminary results suggest however that the 
difference between frozen and cutoff regularisations is of the same order as the width of the inner band. 
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that it is not incompatible with the available structure function data. 

At very small- x one finally sees the BFKL growth of the splitting function. We recall 
that the branch cut, present for a fixed coupling, gets broken up into a string of poles, with the 
rightmost pole located at u c , to the left of the original branch point (u s ), lo s — uj c ~ b^aT [10j. 
The origin of this correction is similar to that of the b 2 ^ 3 a^ 3 contributions to up for cutoff 
regularisations |3Z|. The dependence of u c on Q is shown in Fig. H2 together with its scale 
and IR regularisation dependence. It is slightly lower than the earlier determination in the 
cu-expansion JU| 8 . Both determinations are substantially below ou s , as expected. 

4 Conclusions 

In this letter we have outlined an approach to renormalisation-group improved NLL BFKL re- 
summation that is convenient for numerical determination of the high-energy Green's functions 
and splitting functions. This is an important step on the way to complete RGI NLL predictions 
for high-energy cross sections and small- x structure function evolution. 

We have discussed several important new results obtained within this approach. Most 
striking is the increase in the domain in k and Y that becomes perturbatively accessible once 
one includes subleading corrections. Even the amount of growth (number of 'e-folds') of the 
Green's function that can be calculated perturbatively increases significantly. 

Another result concerns the size of preasymptotic effects: for example, for the transverse 
scale studied in Fig. BFKL growth sets in only for rapidities greater than 4. Hence it is vital 
to study the full Green's function rather than just the high-energy exponents discussed so far 
in the literature, and to include the physical impact factors (THUH] along the lines suggested 
in PS]. 

In the collinear region, a key result is the practical demonstration of factorisation of the 
small-x Green function, and the extraction of the resummed P gg splitting function. Here 
we obtain the high-energy exponent u> c , but we find that preasymptotic effects are again of 
fundamental phenomenological importance. As has been observed in an approach without 
renormalisation-group improvement |13j . the main feature in today's accessible x range is a 
small-a; dip rather than growth (our dip is however much shallower). This phenomenon has yet 
to be fully understood, because several subleading effects are likely to come into play. 

On the whole, the present work encourages us to trust resummed perturbative calculations 
for next generation accelerators, and shows that subleading contributions not only decrease 
high-energy exponents, but provide significant preasymptotic effects. 
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